Giant Skyrmions Stabilized by Dipole-Dipole Interactions 
in Thin Ferromagnetic Films 
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Motivated by a recent magnetization reversal experiment on a TbFeCo thin film, we study a topological 
excitation in the anisotropic nonlinear sigma model together with the Zeeman and magnetic dipole-dipole inter- 
actions. Dipole-dipole interactions turn a ferromagnet into a frustrated spin system, which allows a nontrivial 
spin texture such as a giant skyrmion. We derive an analytic formula for the skyrmion radius. The radius is 
controllable by the external magnetic field. It is intriguing that a skyrmion may have already been observed as 
a magnetic domain. A salient feature is that a single skyrmion can be created or destroyed experimentally. An 
analysis is made also on skyrmions in chiral magnets. 



Skyrmions, originally proposed to account for baryons in 
nuclear physics, now play crucial roles in almost all branches 
of physicslll]. In condensed matter physics, in particular, the 
notion of skyrmion has proved to be most successful, mani- 
festing itself in various experimental observations. Skyrmions 
are solitons in a nonlinear field theory characterized by the 
topological quantum number. In addition to its topological 
stability, for a soliton to materialize, the Hamiltonian must 
contain extra terms that introduce scale to the system. Indeed, 
Skyrme introduced a term of the fourth order in derivatives 
into his original model||2[]. In an instance of quantum Hall 
ferromagnets, the scale of a skyrmion is determined by the 
competitive interplay between the Zeeman and the Coulomb 
interactions yl |4J]. Recently, a skyrmion lattice[5[] as well 
as a single skyrmion[6] have been observed in chiral mag- 
nets, where the role of the Dzyaloshinskii-Moriya interaction 
(DMI) has been discussedJ3,lll. 

In this paper we propose a new mechanism of skyrmion ma- 
terialization based on the magnetic dipole-dipole interaction 
(DDI). Magnetic DDIs exist in all magnetic materials. They 
introduce frustration into a ferromagnetic state. Frustration 
between DDIs and exchange interactions could lead to rich 
phenomena in magnetic materials. Indeed, the ground state of 
a thin ferromagnet film has a complicated magnetic domain 
structure comprised of stripes, bubbles or labyrinths M9I- 1 1311 . 
We show that a giant skyrmion emerges in external magnetic 
field. The size of the skyrmion is so large (~ 1/ito) that it may 
be observed simply as a magnetic domain. 



We present a skyrmion spin texture in the anisotropic non- 
linear sigma model. Deriving an analytic formula for the ra- 
dius of a giant skyrmion, we find that there exists a minimum 
radius. As the magnetic field increases, the radius decreases 
and suddenly shrinks to zero, where a skyrmion disappears. 
The skyrmion spin texture would make a quantum jump to 
the homogeneous ground state, which breaks the conserva- 
tion of the topological quantum number. Furthermore, a sin- 
gle skyrmion can be created by destroying the magnetic order 
within a tiny spot, as we soon argue. It is also pointed out that 
the skyrmion excitation energy becomes negative, making a 
formation of a skyrmion lattice possible, when the external 
field is less than a certain critical value. A skyrmion enjoys 
the topological stability, as far as the external magnetic field 



FIG. 1. (Color online) (a) Illustration of a giant skyrmion (~ 1/im) 
in a thin ferromagnetic film. The simplest spin texture has naturally 
a nontrivial Pontryagin number. It can be created by applying fem- 
tosecond optical pulse irradiation focused on a micrometer spot and 
thus destroying the magnetic order locally, (b) Illustration of mag- 
netic flux lines around a skyrmion due to magnetic dipoles. When the 
magnetic order is destroyed locally, a new order is generated which 
is opposite to that of the environs, so that the magnetic flux closes by 
itself as short as possible. 



is not too large or not too weak. 

A giant skyrmion may have already been observed exper- 
imentally as a magnetic domain in a TbFeCo thin film lll4fl . 
In this experiment the ground state is a homogeneous spin- 
polarized state in a small external magnetic field. By applying 
femtosecond optical pulse irradiation focused on a microme- 
ter spot, it is possible to destroy the magnetic order locally. 
Then, the DDI generates an effective magnetic field, leading 
to a new magnetic order with reversed magnetization inside 
the spot. The resultant spin texture must be that of a sin- 
gle skyrmion, since it is the simplest and most natural spin 
texture, as illustrated in Figfl] Thus, a skyrmion can be cre- 
ated at any point from the ground state experimentally by de- 
stroying the magnetic order within a tiny spot. This process 
breaks down the continuity of the classical fields, and the con- 
servation of the topological quantum number is lost. A giant 
skyrmion may also occur in chiral magnets, where the spin 
twists around the skyrmion so as to break the chiral symme- 
try, as has been illustrated in Refs.|5|-|7|]. 

Our system is a two-dimensional ferromagnetic plane in 
perpendicular magnetic field. The basic Hamiltonian is the 
nonlinear 0(3) sigma model Hj with easy axis anisotropy 



Hj = ~T / d 2 x[(d k n) (d k n) - r 2 (n z )% (1) 



where T = (l/2)zS 2 J is the exchange energy! 10] (z denotes 
the number of nearest neighbors, S the spin per atom, J the 
exchange constant), £ is the single-ion anisotropy constant, 
and n = (n x ,n y , n z ) is a classical field of unit length. The 
ground-state solutions of Hj are n = (0, 0, ±1). 
The DDI between two spins is described by the term 



H 



D 



4tt 



d 2 xd 2 x' 



n(x) ■ n(x') 3[n(x) ■ r][n(x') ■ r] 



(2) 

with r = x — x' and r = \r\, where il = NS 2 g 2 fi B /j,o/a i 
is the DDI strength! 10] (N denotes the number of the layers, 
g the Lande factor, [is the Bohr magneton, a the lattice con- 
stant). 

We apply the magnetic field h perpendicular to the plane, 



H z = -A- 



f d 2 x . . 

/ —j-n z [x), (3) 



with the Zeeman energy Az = Sg/iB^oh, so that the ground 
state is the spin-polarized homogeneous state, n = (0, 0, 1). 

The ordinary ferromagnet is described by the Hamiltonian 
H = Hj + Ho + Hz- On the other hand, the chiral magnet is 
described by the Hamiltonian H = Hj + Hd + Hz + Hum, 
where 



Hum = Dn(x) • (V x n{x)) 



(4) 



is the DMI term. It breaks the chiral symmetry explicitly. 

The use of a continuum approximation and of classical 
fields to represent the spins is justified as far as we analyze 
phenomena whose characteristic wavelength is much larger 
than the lattice constant. 

We first analyze the ordinary ferromagnet and then the chi- 
ral magnet. We start with the study of the ground state in 
weak external magnetic field. It is well known that the DDI 
makes the ground state to have an alternating up-down stripe- 
domain structure. Our results are consistent with those in 
literature! 9- 13] but we a ls° present some new formulas. The 
basic object in our analysis is the one-dimensional solution of 
the anisotropic nonlinear sigma model ([TJ independent of the 
y coordinate. It has two types of solutions; a kink localized in 
the x-axis and a soliton lattice periodic in the x-axis. 

The kink solution of the nonlinear sigma model (Q]i is given 
by n x (x) = \J\ — a 2 {x), n y — 0, n z (x) = a(x) with 



<j(x) = ±tanh[(.T — Xo)/£]. 



(5) 



It describes a domain wall with thickness £ lying along the 
y-axis and located at x = xq, and separates the two field con- 
figurations n = (0, 0, 1) and n = (0, 0, —1). The exchange 
energy is E^ nk — 2TL y /^, where L y is the length of the kink 
along the y-axis. The DDI energy (O is found to be negative 
and divergent for a fixed value of L y , indicating the ground 
state is a condensed phase of domain walls in the absence of 
the Zeeman effect. 

Such a phase is described by the periodic solution, 



<t{x) = sn (#/«£, k 2 ) , 



(6) 




FIG. 2. (Color online) (a) Illustration of a soliton-lattice solution 
a(x) described by the Jacobian elliptic function, (b) Illustration 
of the ground-state structure with an alternating up-down stripe- 
domains separated by domain walls, (c) Illustration of a skyrmion 
field configuration cr(x) described by d 1 51 ) with (16). (d) Illustration 
of a skyrmion spin texture with up and down domains separated by a 
circular domain wall with thickness £. 



in terms of the Jacobian elliptic function with < k 2 < 1: 
See Figflja). The periodicity of a(x) is 2£$ with 



i s = 2 K £K(k 2 ), 



(7) 



where K(n 2 ) is the complete elliptic integral of the first kind. 
The soliton-lattice solution (0 implies that the ground state 
has an alternating up-down stripe-domain structure separated 
by domain walls, as illustrated in FigfSfb). The width of one 
stripe is given by £$■ 

Substituting © into the Hamiltonian (JTJ we obtain the en- 
ergy of the soliton lattice to be 



E L 



2£ 



2 LxLy 



2 E(k 2 ) 
~^ K{k 2 ) 



(L x /£s)E ] 



kink' 



(8) 

where E(k 2 ) is the complete elliptic integral of the second 
kind. The number of solitons is given by L x /£$- The last 
equality holds in the regime where k « 1 or £s/£, ^ 1- 
We are able to estimate the DDI energy for £s/£, 3> 1 as 



^sl = - - 



' ' L'x'-'y 



hi 



is 



2tt £s edf 



(9) 



where dp is the thickness of the film. The Zeeman interaction 
does not make any contribution. 
The total energy is 



E : 



SL 



L X Ly 

~ls~ 



2r n £ s 

£ 2-k edF 



Minimizing this with respect to £s, we obtain 

£ s = d F cxp (2 + 47rl7£ft) . 
The ground-state energy is 

E Sh = -{L x L y n/2Trd F e 2 )e~^ r ^. 



(10) 



(11) 



(12) 



The stripe-domain structure is not robust and may be de- 
formed into a labyrinth structure. The formula ( fTTT i gives the 
order of the width of a stripe segment in this case . 



On the other hand, in strong external magnetic field we ex- 
pect the homogeneous state, which has only the Zeeman en- 
ergy, -Ehorao = — {L x L y /a 2 } Az- By comparing this with 
(fl2l l. there exists the critical magnetic field h c with the asso- 
ciated Zeeman energy being 



( a )tE, 



A| = (a 2 n/2ird F e 2 )e- 4 * r ^. 



(13) 




The ferromagnet phase appears for Az > A z , or h > h c . 

We proceed to analyze a skyrmion spin texture in the fer- 
romagnet phase. In the regime where the classical-field ap- 
proximation is valid, there exists the topologically conserved 
charge, that is the Pontryagin number, 

Qsk y = -^5I j d 2 xe ij n{x)-{d l n{x) x d 3 n{x)) , (14) 



where i,j run over x, y with e,j being the completely anti- 
symmetric tensor. A spin texture possessing nonzero Q s k y is a 
skyrmion by definition. 

We consider the spin texture given by 



h c 2e" 1 /l s h s 



FIG. 3. (Color online) (a) The skyrmion size R is determined by min- 
imizing the energy E s ^ y (R). Various curves are for different mag- 
netic field h. Skyrmion is stable only for h < h s . (b) The skyrmion 
size R decreases as the field h increases. A skyrmion lattice may 
appear for h c < h < 2e~ 1 h s . 



The skyrmion is not physical as it stands, though it carries 
the topological number, since the energy can be made zero 
by letting R — >• in ( fTSI l. It is made physical by the DDI in 
addition to the Zeeman effect, as we now show. 

The DDI energy is estimated by substituting the skyrmion 
configuration ( fToT ) into the DDI Hamiltonian ©, 



= -^l~a 2 (r)cos(9 + 9 ), 



Eg y = -n[Rln(R/d F ) - R]. 



(19) 



-y 1 — a 2 (r) sin (9 + 9$) , n z = cr(r) (15) The total excitation energy is given by 



in the cylindrical coordinate, where 9$ is a constant parameter 
representing a zero-energy mode in the nonlinear sigma model 
Hj, while a(r) is a trial function. The boundary conditions 
are a(r) — > 1 as r — > oo, and a(r) = — 1 at r = 0. The first 
condition implies this describes a soliton in the ground state 
with a(r) = 1. The second condition is required to remove 
the multivalueness of the field n(x) at the center. Substituting 
( fTBI l into (O, we find that Q s ^ y = 1 for any 0q, where the use 
of the boundary conditions on a(r) is made. Hence, the spin 
texture ( fTBI ) has one unit charge, and hence it is a skyrmion. 
Its stability is guaranteed topologically. 

The chiral symmetric configuration is given by the choice 
of 9q — 0, which minimizes the DDI energy. This is the 
simplest and most natural field configuration in the nonchiral 
magnet, as illustrated in Fig[T] 

We consider a domain with radius R separated by a domain 
wall, where a{r) f» — 1 for r < R and cr(r) « 1 for r > R. 
The wall must have thickness £ due to the anisotropy term in 
([T). Such a domain structure is well described by [Fig|2jc)] 



£k y (i?) = ^^ - n[Rln ^-R]+ ^A z , (20) 
£ dp Or 

where the last term is the Zeeman energy. We have illustrated 
Esky(R) for typical values of Az in Figf3] The variation of 
(EOi l with respect to R yields, 



47rr „, R 2ttA z 

—--n\n— + ^-R = 0. 

£ dp a z 



Solving this for R we find 



R=- 



2vrA. 



-W-i 



A; 



(21) 



(22) 



where A| = e{a 2 Q</2iid F e 2 ) exp {-AirT/Qg), and W-x(z) 
is the Lambert function. We may also solve d2Tb for fi, 



n 



47rr/£ + 2irA z R/a 2 
log(R/d F ) 



(23) 



a(r) = tanh [(R/g) \og(r/R)} 



(16) 



which satisfies the boundary conditions at r = and r — >• oo. 
Near the domain wall, it behaves as 



a{r) = tanh[(r - R)/g\ at r « R, 



(17) 



This is interpreted as the formula to determine the effective 
DDI strength from the size of a skyrmion. A comment is in 
order: It follows that A s z = eA c z with ( fl~3"l >. This is a relation 
peculiar to the nonchiral magnet. 

The Lambert function W-i(z) has a real value only for z > 
— 1/e, as implies 



provided R 3> £, as agrees with the domain-wall solution (O. 
Provided R 3> £, a(r) « ±1 except for the domain- wall 
region, and the exchange energy is estimated as 



A z < A|, or h < h s 



(24) 



Namely, the skyrmion is stable only when the magnetic field 
is less than the critical one h s . The minimum radius is 



Ei y = AnTR/C 



(18) 



Rn 



a 2 n/(2TrA z )=ls/e, 



(25) 



where £$ is the width of one stripe in the absence of the mag- 
netic field and given by ( fTTT ). 

When h exceeds h s , the radius shrinks to zero and the 
skyrmion disappears, which is beyond the classical-field 
regime of the spin. What actually happens would be that the 
spin direction at the skyrmion center makes a quantum jump 
to that of the homogeneous ground state. 

The skyrmion energy is less than the ground-state energy 
for h < 2e~ 1 h s , as illustrated in Figf3] Skyrmions would 
condense and form a lattice to make a new ground state when 
the sample is cooled from high temperature in the magnetic 
field h < 2e~ 1 h s . On the other hand, once a skyrmion is cre- 
ated in the ferromagnetic ground state (h > 2e~ 1 h s ) and then 
the magnetic field is decreased at low temperature, it remains 
to be a stable object even for h < 2e~ 1 h s . 

We next discuss a skyrmion in chiral magnets such as MnSi 
and FeCoSi thin films. The Hamiltonian is given by H = 
Hj + H D + H Z + H DM , where H DM is the DMI term ©. The 
system is in the ferromagnetic phase beyond a certain critical 
magnetic field h c . The skyrmion spin texture is given by ( TT5l ). 
where the zero-energy mode Q$ of the nonlinear sigma model 
Hj is fixed by minimizing the total energy. We find 6q = 
— 7r/2 due to the DMI term. Namely, the spin twists around 
a skyrmion as in the illustrations given in Refs.yHZfl. All the 
rest analysis can be carried out almost as it was. Indeed, the 
function a{r) is given by ( [Tol l. Provided R 3> £, the DMI 
energy is calculated to be 



£™ = -2-kDR. 



(26) 



There exists a small increase of the DDI energy due to the 
twisting of spins inside the domain wall, which we neglect. 
As a result the net effect is to renormalize the exchange energy 

as 



r^rv 



r 



1 



-£D 



(27) 



in the skyrmion excitation energy (1201 ), Then, we can derive 
d2l1i^(f25b just as they are simply by replacing T with Tdm- 
A giant skyrmion stabilized by the DDI may also occur in the 
chiral magnet provided T > ^£,D. The essential role of the 
DMI is to twist the spin texture around a skyrmion. 

The present work was motivated by a recent magnetiza- 
tion reversal experiment! 14] on a TbFeCo thin film. It is 
not certain if our analysis can be applicable to this material 
as it stands, because TbFeCo is amorphous and what is real- 
ized is a ferrimagnet. Nevertheless, it is tempting to make a 
comparison of our analytical results to the experimental re- 



sults by identifying the observed magnetization reversal do- 
main as a skyrmion. Sample parameters ! l4l 1511 are the lattice 
constant a = 0.3 [nm], the layer thickness dp = 35 [nm] 
and the exchange energy T — 3.5 x 10 _20 [J]. The experi- 
ment was carried out in the magnetic field at 170[Oe], where 
the radius of the observed skyrmion is about R = 1 [/mi] 
and the thickness of the domain wall is about £ = 25[nm]. 
We then use d23l to determines the effective DDI strength, 
£1 = 8.4 x 10~ 12 [J/m]. This is a good agreement with 
Q, = 7.0 xlCT 12 [J/m] for TbFeCo, which is evaluated from its 
definition given below (|2). We predict the minimum skyrmion 
radius to be R m { n — 0.7 [/an], Using these parameters, the 
stripe width becomes of the order of £$ = 1.9 [/mi] in the 
absence of the magnetic field. On the other hand, h s = 190 
[Oe], beyond which a skyrmion is expected to disappear. 
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